£'-band excitations in the magnetic Keplerate molecule Fcso 
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The low-temperature excitations in the magnetic Keplerate molecule Feso as calculated by linear 
spin- wave theory (SWT), modified linear SWT, and spin-level mean-field theory (SLMFT), are 
compared to the recent inelastic neutron scattering results by Garlea et al. [Phys. Rev. B 73, 
024414 (2006)]. SLMFT reproduces a part of the experimental spectrum rather well, but not all of 
it. SWTs yield a small fraction of the i5-band excitations and hence are not capable of a complete 
description of the excitation spectrum. 

PACS numbers: 33.15.Kr, 71.70.-d, 75.10.Jm 



The subject of the elementary spin excitations in fi- 
nite antiferromagnetic (AFM) Heisenberg spin clusters 
has become of much interest recently, as for many exper- 
imentally available molecular nanomagnets, such as the 
AFM wheels or the Keplerate molecule Feao, the Heisen- 
berg spin Hamiltonian 



H 



(1) 



2J 



is the appropriate starting point for a discussion.^'^ Jy 
measures the exchange interaction between spins i and j, 
and Si is the length of spin i. The number of spin centers 
in the cluster will be denoted as TV. 

For a number of AFM Heisenberg clusters, the low- 
lying energy spectrum exhibits a rather particular struc- 
ture: As function of total spin S, the lowest-lying levels 
form a set of rotational bands (RBs) for which E{S) cx 
S{S + The lowest set of bands is denoted as L 

band; the higher-lying ones as E band (and the remain- 
ing states as quasi continuum) The L band is related 
to a quantized rotation of the Neel vector. For bipar- 
tite systems, the E band corresponds to (discrete) AFM 
spin-wave excitations. This picture of the excitations has 
been confirmed experimentally for a number of bipartite 
systems^, in particular the AFM wheelsji 

Also for the molecule Fcao a RB structure was 
conjectured This system consists of an icosidodecahe- 
dral arrangement of 30 spin-5 /2 centers and is character- 
ized by three AFM sublattices.^ For the existence of the L 
band in this molecule solid evidence from both theory and 
experiment is available Concerning the E band, how- 
ever, the situation is less clear. This brief report is moti- 
vated by two recent studies, a study of the excitations us- 
ing modified spin- wave theory,-^ and an inelastic neutron 
scattering (INS) experiment, where the observed scatter- 
ing intensity was related to the i?-band excitations^^ In 
the following, the spin excitations in Feao as calculated 
by linear spin- wave theory (LSWT), modified linear spin- 
wave theory (mLSWT), and spin level mean-field theory 
(SLMFT), as well as observed experimentally by INS are 
compared. 

For the SWT calculations, Ref. [lli was closely fol- 
lowed, albeit with a different numerical implementa- 



tion. The classical ground state of H (Ref. [Tsh is de- 
termined by the conditions x ^ Jij^j = for each 



which are solved numerically by iterating = 
Si X)j \ X)j until convergence Local coor- 

dinate frames are introduced at each spin site such that 
the local z-axes coincide with the directions of the clas- 
sical spins. Writing Si = {S^ , S~ , Sf )'^ for the spherical 
spin operators in the local frames, the Hamiltonian be- 
comes 
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Expressions for the elements U^J^ {v,^ = -\-,—,z) are 
given in . The spin operators (in the new frames) are 
expressed by Holstein-Primakoff bosons, w \/2Siai, 
S- = Si - ajai, yielding Hb = - Y^i-, JijC^h with 



C. 



(2) 



{S,S, - S^aja, - S^a]aj)U.^f + 2^/S^ 



This real-space bosonic Hamiltonian is Bogliubov diag- 
onalized numerically. We follow Ref. [l3| and introduce 
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and the NxN matrices A and B, which depend on Si, 
Jij and U^J^. The properties = A and B"^ = B ensure 

the hermiticy of D, and Hb- The Bogliubov diagonal- 
ization (or, in the language of p7j . para-diagonalization) 
introduces N new bosons Ck via c^^ ~ a^V^. The 
2Nx2N matrix V obeys V^'iV = i with 1 = diag(l,-l) 
(1 is the NxN unit matrix). This ensures the boson 
character of the Ck- Para-diagonalization corresponds 
to finding a V with (V^)"! • D • V"^ = iE, where 
E = diag(£;i,£;2, . . . ,£;Ar,£;i,£;2, . . . ,£'Ar). The para- 
diagonalized Hamiltonian then reads 
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Numerical algorithms for determining E and V are given 
in The procedure so far corresponds to linear SWT 
(LSWT), i.e., H^j^^^^^ = Hg, The ground-state and 
excitation energies are given by E^'^^^'^^ = -Eo+j J2k 
and E^f'^^'^^ = Ek, respectively. 

Standard SWTs, such as LSWT, start from the as- 
sumption of an ordered, symmetry-broken ground state, 
which is obviously incorrect for finite clusters. In the 
modified SWTs,^^'^^'^^ spin- rotational invariance is re- 
stored "by hand" by enforcing zero on-site magnetiza- 
tions, (Sf) = Si — {a\ai) = ((.) denotes the quantum 
expectation value in the ground state) . N Lagrange mul- 
tipliers are introduced, and the boson Hamiltonian of 
modified hnear SWT (mLSWT) becomes 
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Hb + ^^^ii{a\a., - Si). 



(6) 



With A', 
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{fii/2)5ij and E'q 
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assumes the form of 



J2iMSi + 5), also H}^ 
Eq. ([3]). Para-diagonalization yields H 
J2k E'kic^c'f. + \)- The ground-state and excitation ener- 
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gies, 4™^^^^) ^E', + \Y.kE'k and i?^'";-" ^ > = 
implicitly depend on the Langrange multipliers, which 
are found numerically by solving jj'^^^^^^) iteratively 
until {a\ai) = Si for each i. 

A further approach, called spin-level mean-field theory 
(SLMFT), was explored. The RB structure assumed for 
Peso implies approximating the wavefunctions of the L 
band by the spin levels \aSASBScSM) , where 5*^, Sb, 
and Sc are the total spins of the three AFM sublattices, 
and Sa = Sb = Sc = SiN/3 {a denotes intermediate 
spin quantum numbers) i^i^i^i^S This is equivalent to de- 
scribing the L band by the three-sublattice Hamiltonian 
Habc = -6J/N{Sa -Sb + Sb-Sc + Sa- Sc). 

The E band, which corresponds to excitations out of 
the L band, is then related to the spin levels with Sa 
— Sb — SiN/3, Sc = SiN/3—1, and the permutations 
thereof i^i^i^S These states are degenerate for Habc, but 
are expected to be split by the perturbation H — Habc- 
This suggests, as a first approximation, to diagonalize 
H in the subspace of the spin levels IuSaSbScSM), 
with 5*^, Sb, and Sc chosen as appropriate for the 
L and E bands. In this work, where only the zero- 
temperature excitation spectrum is considered, H was 
diagonalized in the subspace of the S—0 ground state 
and the S=l states of the L and E bandsi^i Accord- 
ing to the coupling rules of three spins, the S=0 ground 
state is non-degenerate, and the 5*=! sectors of the L 
and E band consist of 3 and 81 spin levels, respec- 
tively. The numerical calculation of the matrix elements 
{aSASBScSM\H\aSASBScSM) is simple thanks to the 
irreducible tensor operator techniquesi^ 

Sublattice Hamiltonians, such as Habc, may be "de- 
rived" by Fourier analysis of the full Hamiltonian^iii^ 
This procedure works for spin clusters with large sym- 
metry, such as rings, polyhedrons, or extended lattices. 
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FIG. 1: (a) Energy levels of Feso as calculated by LSWT, mL- 
SWT, and SLMFT (gs = ground state; numbers in brackets 
give the degeneracy), (b) T — neutron scattering intensity 
calculated by SLMFT. The intensity of the peak at 0.2| J| was 
divided by a factor of 40. 



It may be constructed also by replacing each spin by the 
mean-field spin = J2ier] corresponding 
sublattice ry, which provides a more general method. ^'^ 

The results of the calculations for the excitations in 
Feso, as obtained by LSWT, mLSWT, and SLMFT, are 
displayed in F ig, [p a) (the mLSWT results were obtained 
previously in [il|). In LSWT, the excitation spectrum 
consists of 3 levels at zero energy (with respect to the 
ground state), followed by 8 levels in the energy range 
0.0078 to 0.0247| J|, and four bands in the energy range 
3.5— 5.5| J|. As expected, mLSWT strongly affects the 11 
lowest excitations, the spectrum is in particular gapped. 
The first 10 excitations are at 0.2758| J|, followed by one 
at 0.3900| J|. The higher-lying excitations are affected by 
less than 0.3% (the levels at ca. 3.5 |J| are weakly split 
into two subgroups). Clearly, the 11 low-lying excita- 
tions, and the ground-state, should be related to the L 
band, while the remaining higher-lying excitations should 
be related to the E band. SLMFT produces a more struc- 
tured energy spectrum. It consists of the three triplets of 
the L band at 0.2| J|, and the 81 triplets of the E band, 
which span the range 3.7811 to 7.3682|J|. For compari- 
son, Habc gives energies of 0.2 and 5.2| J| for the L and 
E band, respectively. 

All three theories yield a significant splitting of the 
E band, which is about 2\J\ for the SWTs and 3.5|J| 
for the SLMFT. The experimental INS data at 65 mK, 
which is reproduced in Fig. [21 in fact shows a broad fea- 
ture from 0.2 to 1.1 meV. Within the three-sublattice 
Hamiltonian approach it was modeled by a single Gaus- 
sian of width 0.66 meV.^^ Its position was determined 
to 0.56 meV, which corresponds to J = —0.108 meV, 
in rough agreement with the value J = —0.134 meV in- 
ferred from magnetic susceptibilityjiii^ The above the- 
ories partially account for the observed broadening (as 
noted before for the SWTs in [Oll^). Since SLMFT pro- 
duces the largest splitting, it is compared to experiment 
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FIG. 2: Inelastic neutron scattering spectrum for Feao. The 
solid circles show the T = 65 mK experimental data of [l^ 
(data were recorded at OSIRIS with final neutron energy Ef 
= 1.845 meV and integrated over the range Q = 0.9—1.8 A~^). 
The solid line and the vertical lines represent the T = scat- 
tering intensity as calculated from SLMFT (J = —0.125 meV; 
the solid line was obtained by a convolution with a Gaussian 
of linewidth 0.27 meV). The inset displays the experimental 
data after substraction of the theoretical curve. The arrow 
marks intensity not accounted for by the theory. 



in more detail. The powder INS intensity integrated over 
Q = 0.9-1.8 was calculated from the wavefunctions 



obtained by SLMFT using the formulas of [2j| and the 
crystallographic Fe positions of Feso-— The result for J = 
—0.125 meV is shown as the vertical lines in Fig. [2l De- 
convolution of this spectrum with a Gaussian of width 
0.27 meV produced the solid curve shown in Fig. [2l The 
theoretical spectrum exhibits a strong peak at 0.025 meV 
due to the excitations from the ground state to the L- 
band triplets. It is buried under the large quasielastic 
contribution and hence not detected in the experiment. 
Interestingly, the experimental data in the range 0.4 to 
1.1 meV is very well reproduced by SLMFT, and the ob- 
tained value for J agrees within 7% with the value from 
magnetic susceptibility, which is a significant improve- 
ment over the value obtained with Habc- However, in 
the range 0.2 to 0.4 meV there is additional intensity. It 
might be described by an "appropriate" background, but 
considering the resolution function of OSIRIS it should 
be assigned to a feature not correctly reproduced by the 
SLMFT. In fact, subtracting the theoretical curve from 
the experimental data (see inset of Fig. [2]) strongly sug- 
gests additional scattering from Feao in this range. 

For the discussion, first the conceptional differences be- 
tween the SWTs and SLMFT are noted. The SWTs start 
from a symmetry-broken ground state, while SLMFT 
works with spin levels and is intrinsically spin rotational 
invariant. Hence, the relation between the states ob- 
tained by SWTs on the one side and SLMFT (or exact 
results) on the other side is not always obvious. For a bi- 
partite system, this is not a big problem. SWTs produce 
N excitations which relate to N— 1 spin excitations in the 
5=1 sector, i.e., one L-band and A^— 2 £'-band triplets 



(for a system with S=0 ground state). There is hence an 
obvious correspondence. The L band is associated to the 
(quantized) rotation of the Neel vector, and the E band 
to the (discrete) spin-wave excitations j^i^i^iii^ 

For a system with three (or more) sublattices the sit- 
uation is less clear. The SWTs again produce N exci- 
tations, while SLMFT concerns with 3 L- and 3(7V— 3) 
i?-band triplets (for a system with three sublattices and a 
5=0 ground state). Clearly, since there are many more of 
them, not all of the £'-band triplets can be spin-wave ex- 
citations. This has an important implication. The (zero- 
temperature) spin-correlation functions are governed by 
the transition matrix elements (0|S'j'^|n) (|0) denotes the 
ground state, |n) the nth excited level). Within the 
sublattice-Hamiltonian approximation therefore only the 
transitions from the ground state to the triplets of the 
L and E bands have non-zero intensity, while all other 
transitions have zero intensity!^ That is, only the triplets 
of the L and E bands are relevant. On the other hand, 
there is no general restriction on the intensities of these 
transitions. That is, all triplets of the E band are rele- 
vant. This is underlined by Fig. [ijb), which reproduces 
the INS intensity as obtained from SLMFT. Comparison 
with Fig. [T{a) shows that all parts of the energy spec- 
trum contribute. Since the number of E-hand triplets 
considerably exceeds the number of magnons, one nec- 
essarily has to conclude that (for three and more sub- 
lattices) SWTs are intrinsically not capable of describing 
the spin dynamics correctly, or completely, respectively. 

These points are in favor of SLMFT. However, a com- 
parison of SWTs, SLMFT, and exact results for even- 
numbered AF rings revealed that SLMFT provides an 
overall estimate for the spin-wave energies, but fails to 
reproduce their dispersion (which SWTs of course do)i^ 
On the other hand, the splitting of the E band in AF 
finite chains due to the open boundary conditions is well 
graspedi^l Apparently, SLMFT accounts for some but 
not all features of the energy spectrum. 

Concerning SLMFT, the studies on the rings and 
chains showed the following trends. The upper edge 
of the exact .E-band spectrum was obtained reasonably 
well; it was overestimated by about 10%. The disper- 
sion of the magnons, however, was greatly underesti- 
mated, yielding very narrow E bands: The exact i?-band 
spectrum extends to much lower energies as obtained 
by SLMFT. These trends could explain the discrepan- 
cies in the SLMFT analysis of the Fe^Q INS data. First, 
a J value which is 7% too small is consistent with the 
the trend to overestimate. Second, the scattering inten- 
sity in the range 0.2 to 0.4 meV could be associated to 
magnons, which are at lower energies than calculated be- 
cause SLMFT does not account for their dispersion. 

Concerning the SWTs on the other hand, they were 
found to reproduce the structure of the E band very well 
for the rings and chains, but to underestimate energies, 
i.e., to provide lower bounds for the i?-band spectrum 
(the accuracy of LSWT and mLSWT was ca. -20%, 
that of interacting SWT and modified interacting SWT 
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about —10%). With respect to Feao, Figure [TJ a) then 
would imply that there are no excitations below 3.5 |J| 
(besides the L-band excitations at very low energies), 
contradicting an association of the 0.2—0.4 meV inten- 
sity to magnons. For the triangular AFM lattice, how- 
ever, the inclusion of 3rd- and 4th-order boson terms 
was found recently to result in a substantial suppres- 
sion of the magnon energies Since a similar behav- 
ior is expected for Feao, the experimental intensity at 
0.2—0.4 meV should be indeed associated to magnons, 
consistent with the above conclusion from SLMFT. 

The calculated ground-state energies are — 220.802| J|, 
-220.0951 J|, and -195|J| for the LSWT, mLSWT, 
and SLMFT, respectively. A recent DMRG calculation 
yielded -211(2)| Jj.^^ As observed also for the rings and 
chains, LSWT yields a lower ground-state energy than 
mLSWT, while SLMFT is too high. This is linked to the 
fact that SLMFT does not describe the spin-wave disper- 
sion, and hence misses the related quantum fluctuations. 
That LSWT and mLSWT, however, yields ground-state 
energies well below the DMRG result should be consid- 
ered another hint for the problems of the SWTs. 

In conclusion, three theories, LSWT, mLSWT, and 
SLMFT, have been tested for Feao- The SWTs, as they 
do not cover all relevant excitations of the E band, are 
not able to describe the spin dynamics, as measured by 
INS for instance, correctly. SLMFT was successful to 
some extend. It reproduces a significant part of the ex- 
perimental INS spectrum, and yields a J value which is 
in substantially better agreement with the J value from 
magnetic susceptibility than that obtained with Habc- 
However, it does not account for all the experimentally 
observed excitation intensity. 

For bipartite systems the combination of the results of 
SWTs and SLMFT yields a rather good picture of the 
ground state and low-temperature excitations, even in 
a quantitative sense^^ For the three-(and more) sublat- 
tice systems, these techniques could be considered useful 
to some extend for the analysis of experimental data, in 
particular in the absence of better methods. Conception- 
ally, however, they all are unsatisfactory. SWTs, on the 
one side, break the symmetry, which one has then to cor- 
rect for somehow afterwards. Furthermore, they do not 
reproduce the complete E band. SLMFT, on the other 
hand, while overcoming these two issues, is not able to 
yield the dispersion of the magnons correctly. 

In order to further explore the situation for the three- 



(and more) sublattice systems, a comparison of the 
SWTs and SLMFT with exact results would be crucial. 
The spin-1/2 analogue of Feso can be handled by exact 
diagonalizationJ^ The experience for bipartite systems, 
however, indicates a breakdown of the L / E-hand concept 
for spin-1/2; hence the findings might not be transferable 
to systems with larger spins.= A comparison with the tri- 
angular AFM lattice (TL) is thus suggested, as the con- 
cepts relevant here and for Feso are basically identicaL^i!^ 

Interestingly, the A^=9 TL is exactly solved by Habc, 
with proper values for 5*^, Sb, and Sc (a similarity in 
sequence is noted: for bipartite systems Hab is exact 
for the dimer, square, and iV=8 2D square lattice, for 
three-sublattice systems Habc is exact for the trimer, 
octahedron, and A^=9 TL). Hence, the L/i?-band con- 
cept, and all its implications, is exact here. For rings, 
which may be regarded as extensions of the dimer and 
square, the approximation by Hab becomes the better 
the larger Si and the smaller N.^ For TLs, which may be 
regarded as extensions of the trimer and A^=9 TL, one 
may speculate about a similar trend of the accuracy of 
Habc with Si. 

It would be particularly interesting to study the spin- 
5/2, A^=12 TL, because exact diagonalization should be 
possible (though not on the author's computers). The 
spin-1/2 case was deeply investigated previously. "'■^'^ Here 
it is found that the three L-band triplets are slightly split, 
and that the splitting of the L band increases with in- 
creasing S. This should be expected to happen also for 
the spin-5/2 case, and for Fcsq. This effect could be rele- 
vant to explain the observed unusual field dependence of 
the INS intensity)^ with increasing field the thermal pop- 
ulation of the i-band states would change, so that the 
relative contributions of the ii'-band excitations would 
change too. 
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